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ABSTRACT 


The  Optical  Matched  Filter  Image  Correlator  (OMFIC)  can 
identify  targets  and  their  number,  and  provide  information  about 
their  bearing,  positional  coordinates,  range  interval,  and  velocity. 
Thus,  it  provides  a new  concept  for  the  missile  or  projectile 
guidance  application. 

The  optical  memory  size  and  real  time  transducer  update 
requirements  are  two  of  the  most  important  aspects  of  the  OMFIC 
system  when  used  for  guiding  a projectile.  An  analysis  of  these 
factors  shows  that  the  memory  bank  requirements  are  manageable 
and  within  current  technology  capabilities.  In  order  to  cover 
the  range  from  acquisition  (1.5  km)  to  approximately  50  meters 
from  impact  only  397  memory  positions  are  required  for  highly 
discriminating  matched  filters. 

The  transducer  update  rate  can  be  set  at  television  frame 
rates  for  guidance  to  50  meters  from  impact  at  approximately 
M^  “ 0.5.  Under  these  conditions,  the  correlation  peak  should 
remain  at  c ^-3dB. 

Other  aspects  of  the  OMFIC  memory  bank  and  transducer 
requirements  for  missile/projectile  guidance  are  developed  in 
the  report. 
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As  early  as  1971,  a semiactive  guided  projectile  was  gun 
launched  and  functioned  properly  according  to  one  report  (Ref.  1). 

A semiactive  guided  projectile  is  one  in  which  a selected  target 
is  illuminated  with  a laser  designator  by  a forward  observer. 

The  optical  return  picked  up  by  the  projectile  is  focused  upon 
a detector.  This  system  was  mounted  in  an  eight- inch  projectile 
and  utilized  a quadrant  detector  from  which  a position  proportional 
error  signal  for  guidance  was  derived;  the  field  of  view  was 
36  degrees.  Further  work  was  conducted  on  a five- inch  guided 
projectile  with  a 25  degree  field  of  view  and  which  used 
proportional  navigation. 

Additional  developments  improved  the  laser  designated  guided 
projectiles,  and  in  sane  cases,  special  electro-optical  projectile 
developments  were  also  made;  although  these  are  not  considered 
guided  projectiles,  they  do  represent  efforts  to  provide  target 
identification  for  artillery.  One  of  these  is  a 155  mm  artillery 
shell  with  a charge  coupled  device  (CCD)  television  system  contained 
in  it.  At  last  report  (Ref.  2)  it  had  been  successfully  operated 
in  a simulated  flight.  This  effort  represents  the  growing  interest 
of  artillery  forces  to  develop  a "fire  and  forget"  type  of 
scenario,  i.e.,  a system  that  will  continue  to  guide  the  projectile 
without  human  assistance  such  as  the  forward  observer  who  desig- 
nated the  target  with  a laser  beam  in  the  aforementioned  system. 

The  maturity  of  the  guided  projectiles  effort  is  further  reflected 
in  the  recent  issue  of  a specification  for  a projectile  (Ref.  3). 

Several  more  recent  developments  demonstrated  that  optical 
correlators  might  have  potential  application  in  artillery  shells 
(or  for  that  matter,  missiles)  where  in  addition  to  the  "fire  and 
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forget"  feature  one  could  add  automated  target  identification, 
selection,  and  terminal  guidance.  These  developments  are 
documented  in  Refs.  4 and  5. 

A subsequent  correlator  investigation  (Ref.  6)  has  shown  that 
the  parametric  sensitivities  of  a target  (e.g.,  scale,  orientation) 
are  not  so  severe  as  to  make  target  identification  within  a class 
of  targets  difficult.  Rather,  a limited  number  of  views  of  one 
target  would  be  sufficient  to  identify  the  desired  target  among 
similar  but  not  the  same  targets  (at  any  target  orientation). 

This  report  examines  some  of  the  requirements  necessary 
for  an  optical  matched  filter  image  correlator  (OMFIC  ) to  be 
considered  for  both  identification  and  a fire  and  forget  type  of 
guided  projectile.  Considering  the  amount  of  information  that 
the  OMFIC  system,  in  principle,  can  reveal  about  one  or  more 
targets,  it  would  seem  to  be  a logical  candidate  to  be  considered 
for  this  role.  In  the  next  section,  a brief  description  of  the 
OMFIC  system  is  given  after  which  a general  review  of  the  concept 
for  using  the  OMFIC  system  in  this  type  application  is  given. 

There  are  many  critical  questions  which  the  application  of 
an  OMFIC  system  raises  such  as:  What  is  the  memory  bank 
requirement?  What  is  the  update  rate  for  the  transducer?  and 
many  others.  Specific  answers  are  not  always  possible  without  a 
full  specification  for  the  system  but  the  sections  which  follow 
do  present  an  examination  of  the  questions  and,  at  least,  give 
some  bounds  on  the  problem  from  which  answers  may  be  drawn  for  a 
more  definitive,  prescribed  system. 


Optical  Matched  Filter  Image  Correlator,  U.S.  Patent  No.  3,779,492, 
Grumman  Aerospace  Corporation 
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2.  THE  OPTICAL  MATCHED  FILTER  IM4.GE  CORRELATOR 


The  OMFIC  is  a system  in  which  an  image  is  presented  to  a 
matched  filter  of  the  desired  target  and  a comparison  made  to 
determine  whether  a correlation  exists.  This  is  the  process 
common  to  all  optical  correlators  of  the  type  described  by  van 
der  Lugt  (Ref.  7)  and  others.  There  are  two  aspects  to  the  OMFIC 
that  are  different:  Consider  Figure  1.  First,  the  image  is 
presented  (with  natural  light  or  through  some  assist)  to  a 
transducer  (XDR)  where  the  image  is  stored.  The  image  then 
modulates  a collimated  laser  beam  and  the  modulated  beam  is 
presented  to  a transform  lens.  The  second  feature  of  the  OMFIC 
is  that  the  Fourier  transforming  lens  is  a multiple  holographic 
lens  which  takes  the  Fourier  transform  of  the  input  image,  and 
directs  the  transform  in  replication  to  each  of  the  many  matched 
filters  in  the  memory  bank.  Each  of  the  matched  filters  represents 
a different  aspect  of  one  target,  or  different  aspects  of  many 
targets.  For  example,  the  memory  bank  may  contain  all  views 
necessary  to  cover  all  orientations  of  a target  of  one  size, 
or  all  orientations  of  many  sizes.  For  practical  reasons,  the 
matched  filters  are  made  in  conjunction  with  the  holographic 
lens  used.  For  example,  if  the  lens  is  a 10  x 10  array,  then 
there  would  be  at  least  100  matched  filters.  (There  could  be 
more  — this  questions' is  investigated  in  Ref.  6.)  The  matched 
filters  are  prepared  beforehand,  and  thus,  there  is  an  a priori 
knowledge  of  the  character  of  each  filter.  This  is  important 
because  with  this  information,  together  with  knowledge  of  the 
Fourier  transform  properties,  the  OMFIC  system  does  in  principle 
contain  the  following  information: 

• Identification  of  the  target  or  targets 

• Bearing  of  the  target (s) 
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Fig.  1 Elements  of  Basic  OMFIC  System 


• Positional  coordinates  of  each  target 

• Range  (within  an  interval)  to  each  of  the  target (s) 

• Number  of  targets  (and,  if  desired,  the  numbers 
of  several  types  of  targets) 

• If  provided  for,  the  speed  and  direction  of  the 
targets 'headings 

The  remaining,  though  no  less  important  components  of  the 
OMFIC  system,  are  the  inverse  transform  lens  and  the  correlation 
plane  processor.  The  lens  presents  to  the  processor  the  trans- 
form of  the  product  of  the  matched  filter  optical  transfer 
function  (OTF)  and  the  input  image  OTF.  Thus,  the  correlation 
plane  contains  the  input  image  as  modified  by  the  filter.  The 
effect  of  the  transform  and  inverse  transform  operations  is  to 
cause  a spatial  plane  rotation  of  tt  radians.  This  is  a small 
point  but  it  does  caution  one  to  keep  the  image/correlation 
plane  directions  correctly  oriented. 

For  the  present  application  there  are  two  main  questions 
which  are  to  be  examined: 

• As  the  guided  projectile  approaches  a target,  a 
new  view  of  the  target  must  be  recorded  upon  the 
transducer  every  so  often,  and  a test  through  the 
memory  bank  made;  thus,  how  fast  must  the  image 
on  the  transducer  be  updated? 

• The  particular  missile  or  projectile  considered 
will  dictate  the  range  and  target  requirements  to 
be  met;  thus,  how  many  filters  must  the  memory 
bank  storage  contain? 
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We  will  assume  (for  this  analysis)  information  derived  from 
Ref.  6 as  being  representative  of  a typical,  desirable  target. 
Thus,  the  target  of  interest  selected  is  an  M-60  Al  tank.  It 

is  22  feet  long,  12  feet  wide,  and  10  feet  9 inches  high.  It 
was  shown  in  Ref.  6 that  when  the  scale  of  the  image  was  changed 
by  ± 20  percent  in  area,  the  autocorrelation  peak  in  the 
correlation  plane  was  reduced  by  3 dB.  Similarly,  when  the  image 
was  rotated  ± 8°,  the  autocorrelation  peak  was  reduced  3 dB. 

These  are  called  the  parametric  sensitivities  of  the  target  using 
a high  optical  spatial  frequency  matched  filter  and  they  determine 
the  memory  bank  storage  requirements  when  particular  scale 
ranges,  or  orientation  coverages  are  to  be  achieved. 

For  the  projectile,  a 20°  field  of  view  with  an  incident 
angle  of  70°  is  assumed  at  an  acquisition  range  of  1.525  km.  At 
different  times,  projectile  velocities  of  214,  305,  and  400  m/sec 
are  assumed.  The  first  is  considered  a lower  limit  for  any  case 
while  the  other  two  are  considered  nominal  and  extreme  upper 
limit,  respectively.  These  values  are  based  upon  information 
obtained  from  those  experienced  in  the  ordnance  field. 


3.  OMFIC  GUIDED  PROJECTILE  SYSTEM  CONCEPT 


In  concept,  the  scenario  is  shown  in  Figure  2.  A cannon 
launched  guided  projectile  (GP  - e.g. , 155  mm)  reaches  acquisition 
time  at  apogee  and  has  a field  of  view  which  may,  or  may  not 
include  target (s)  of  interest.  These  targets  may  fall  anywhere  in 
the  observable  field. 

At  the  acquisition  range  turn-on,  the  master  control  unit 
(Figure  3)  initiates  a clock  generator  from  which  all  functions 
are  timed.  This  unit  turns  power  on  and  causes  the  various 
sequences  to  start.  These  include  the  transducer  cycling,  detector 
array  scanning,  laser  cycling,  detection,  and  thresholding  controls. 

With  all  units  powered  and  cycling,  the  shutter  opens  and 
the  transducer  is  powered  so  that  an  image  is  stored.  The 
transducer  is  postulated  as  cycling  at  twice  per  second.  After 
the  read  cycle  is  complete,  the  shutter  is  closed  and  the  laser 
activated.  This  enables  the  stored  image  to  spatially  modulate 
the  collimated  laser  beam.  This  image  is  then  Fourier  transformed 
by  the  holographic  transform  lens  which  presents  the  spectrum 
simultaneously  to  all  positions  of  the  matched  filter  memory 
bank.  At  the  latter,  correlation  of  the  image  takes  place  if  a 
tactical  target  of  interest  is  present.  The  inverse  transform 
of  the  product  of  the  matched  filter  OTF  and  image  OTF  is  taken 
and  presented  to  the  detection  array  which  has  been  activated 
and  scanned.  Activation  includes  battery  power  as  well  as 
threshold  setting.  Signals  exceeding  threshold  activate  the 
target  decision  logic  where  it  is  determined  whether  the  target 
is  a primary,  secondary,  or  tertiary  one.  This  decision  is  based 
upon  information  which  has  been  stored  a priori  in  the  matched 
filter  during  fabrication.  At  the  same  time  that  the  target 
decision  process  is  taking  place,  the  positional  coordinates  are 
being  generated. 
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Upon  the  target  position  decision,  the  offset,  or  guidance 
requirement  to  reduce  offset  is  determined.  This  information 
is  then  used  to  generate  guidance  power  signals  which  effect  an 
appropriate  aerodynamic  surface  correction. 

It  can  be  seen  that  one  of  the  key  elements  in  the  use 
of  the  OMFIC  technology  in  a cannon  launched  guided  (CLG) 
projectile  is  the  memory  bank  requirement.  The  latter  must 
consider  all  the  necessary  tactical  positions  of  primary, 
secondary,  and  tertiary  targets  that  are  necessary  in  order  to 
detect  their  presence  within  the  projectile  field  of  view,  range, 
and  other  aspects  affecting  the  input  image.  The  second  key 
element  is  the  rate  at  which  the  transducer  can,  and  must,  be 
updated  with  a new  input  image.  Each  of  these  elements  is  now 
examined. 


4.  IMAGE  SIZE  CHANGE 


One  of  the  most  critical  aspects  of  the  guided  projectile 
scenario,  as  it  affects  the  design  of  an  OMFIC  guidance  system, 
is  the  change  of  image  size  during  the  trajectory  of  the  pro- 
jectile. 

Consider  the  geometry  shown  in  Figure  4,  where  the  projectile 
is  illustrated  at  the  maximum  acquisition  range.  The  field  of 
view  a is  prescribed  at  20°,  and  at  acquisition,  the  projectile  is 
observing  the  scene  from  an  incident  angle  9 ■ 70°  measured 
from  the  surface  normal.  From  the  figure  we  have 

optical  centerline  slant  range  in  km 

outer,  inner  ranges  at  edges  of  field  of  view 
in  the  plane  normal  to  the  ground  surface  and 
containing  missile's  path,  in  km 

a field  of  view,  degrees 

9 incident  angle,  degrees 

F (MFIC  system  input  optic  focal  length,  mm 

imaged  target  (linear)  size,  mm 

S actual  target  size,  mm 

Also,  from  the  geometry  shown  in  Figure  4,  information  in 

L R 

Table  1 might  also  be  obtained  where  X ’ represent  the  perpen- 
dicular distances  off  the  normal  plane  from  the  target  on  the 
gr  ound . 
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TABLE  1 GEOMETRICAL  QUANTITIES  AT  ACQUISITION 


Target  Offset 
Target  Altitude 


“ 1.433  km 
H - 0.265 

X'*'  = 1.525 

x’  = 0.530 

X^’^  = 1.549 


The  target  of  interest  is  imaged  upon  the  transducer  with 
a size  that  depends  upon  the  projectile  imaging  lens  focal 
length,  target  size,  and  slant  range.  Along  the  centerline  slant 
range,  the  linear  target  size  is  given  by 


where 

= centerline  slant  range 

F “ lens  focal  length 

S = actual  target  linear  dimension 

Consistent  dimensions  must  be  used  to  obtain  L^  in  mm.  The  range 
R^  can  be  replaced  by  R^  or  R for  the  minimum,  or  maximum  sizes 
in  the  field  of  view  along  the  edges  of  this  field  of  view. 

As  a reasonable  approach,  a good,  commonly  used  focal  length 
of  305  mm  will  be  generally  employed  except  where  indicated.  Also, 
a nominal  value  for  centerline  slant  range  for  one  projectile 
is  1.525  km  at  acquisition.  These  two  values  will  be  considered 
representative  and  used  in  this  analysis. 
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Thus,  using  a 305  nim  focal  length  imaging  lens,  the  linear 
image  size  for  the  given  slant  range  and  20  degree  field  of  view 
and  depression  angle  yields  the  values  shown  in  Table  2 for  a 
target  such  as  the  M-60  Al  tank  which  is  22  feet  long.  Correction 
has  been  made  for  the  height  (10.75  feet)  of  the  tank.  Also, 
the  scale  factor  referred  to  in  the  usual  aerial  reconnaissance 
sense  is  also  given.  The  change  of  approximately  3 to  1 in  image 
size  between  the  near  slant  range  and  the  far  slant  range  ought 
to  be  noted.  These  values  actually  form  the  limits  in  image 
size  at  acquisition  for  this  field  of  view. 

TABLE  2 IMAGE  SIZES  AND  SCALE  FACTORS  AT  ACQUISITION  RANGES 


slant 

Range 

('km') 

Scale 

Factor 

r" 

1.043 

1.962 

3420 

^o 

1.525 

1.342 

5000 

r'*’ 

2.985 

0.686 

9787 

We  might  note  at  this  point  that  the  initial  or  acquisition 
range  size  of  the  image  impacts  the  required  resolution  of  the 
system;  particularly  the  resolution  of  the  transducer  and  the 
resolvable  detection  element  size  in  the  output  correlation  plane. 
First,  we  consider  the  transducer  requirements.  If  we  were  to 
require  that  the  tank  have  five  cycles  across  its  height  (levels 
determined  in  Ref.  9 for  acceptable  correlative  detection) , then 
for  the  acquisition  range  of  1.5  km  and  a target  height  of  3.279 
meters , we  have 
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5 cycles  x 


rad 

1500M 


s 2.29  cycles/mrad 

If  we  are  to  image  the  entire  field  of  view  of  20°  (or  349 
rarad)  this  would  require  a system  resolution  of  (2.29  cycles/mrad) 
X (349)  = 800  cycles  across  the  processing  aperture.  For  example, 
in  the  transducer  we  must  accommodate  800  cycles  in  the  aperture 
and  if  we  take  a varied  35  mm  aperture  this  requires  a transducer 
resolution  capability  of  23  cycles/mm  which  is  within  the  current 
state  of  the  art. 

The  second  impact  of  the  size  of  the  image  at  acquisition  is 
the  requirement  imposed  on  the  detection  plane;  namely,  that 
ideally  the  detector  must  not  only  respond  to  the  peak  value  of 
the  correlation  signal.  From  a practical  standpoint  we  can 

relax  this  requirement,  in  fact,  based  upon  the  results  from  Ref. 

„ . . correlation  spot  size  . 

9 the  ratio  -rt r : : 1 — — : is  equal  to  0.34. 

object  size  in  correlation  plane  ^ 

Thus,  a corresponding  detection  plane  requirement  calls  for 
0.34  X 800  = 270  lines,  a reasonable  system  in  the  current 
state  of  technology  in  which  about  5 lines  can  be  obtained  across 
the  correlation  spot. 

Returning  now  to  the  scale  change  during  the  trajectory,  we 
note  that  in  carrying  out  our  analysis  we  will  generally  be 
concerned  with  the  three  velocities 


V = 214  m/sec  h Mach  0.64 
305  m/sec  s Mach  0.91 
400  m/sec  s Mach  1.19 
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which 


which  will  be  used  in  conjunction  with  the  ranges  R^,  R ’ which 
are  often  presented  as  the  ratio  (v/R) . Thus,  one  or  more  of 
the  values  of  the  ratio  (v/R)  = 0.14,  0.20,  0.26,  0.32,  and  0.38, 
will  be  used  in  the  analysis  to  bracket  the  combinations  in  the 
ranges  of  interest. 

For  this  analysis  the  path  of  the  projectile  is  straight  to 
the  target.  No  allowance  is  made  for  ballistic  path,  or  guidance 
correction.  During  projectile  flight  there  will  be  a variation 
in  the  size  of  the  image  placed  upon  the  transducer,  and  the 
variation  of  this  image  size,  and  the  frequency  with  which  it  is 
changed  and  held,  is  paramount  to  the  determination  of  the  memory 


bank. 


If  the  projectile  velocity  is  given  by  v,  and  the  acquisition 


slant  range  by  R^,  then  the  image  size  will  vary  such  that,  from 
Eq.  (1) 


n - L ) = ^ ^ 

^ 1 o''  (R  - vt)  R 

' o ' o 


ZS 


O 1 


(f)' 


’image  Size  at'  (v/R  )t 
= Acquisition  : — 

^ / V \ 

.on  TransducerJ  1 - i~lt 


The  image  size  change  factor,  allowing  for  the  initial  input 
image,  then  becomes 


(f) 


(1  - vt) 
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where  V = (v/R^)  is  the  velocity  factor.  Equation  (3)  would 
enable  the  time  history  of  the  image  size  to  be  determined  from 
acquisition  to  impact  (1  - Vt  = 0).  In  Table  3 a computation  is 
summarized  which  shows  how  the  image  change  factor  varies  close 
to  the  impact  time.  The  latter  is  also  given.  In  addition, 
for  each  case  the  image  change  factor  is  given  for  the  running 
time,  50  milliseconds  from  impact.  This  table  should  be  compared 
to  the  curves  in  Figure  5 for  V = 0.14,  0.20  where  the  image 
change  factor  can  be  examined  in  closer  detail  near  acquisition. 
Equivalent  curves  for  the  other  values  of  V would  fall  to  the 
left  of  those  shown,  and  asymptotically  approach  the  respective 
impact  times. 

TABLE  3 IMAGE  CHANGE  FACTOR  NEAR  PROJECTILE  IMPACT 


v/R 

Imp  ac  t 
Time 

Running 

Time 

(sec) 

Image 

Change 

Factor 

X 

0.  14 

7.143 

6.90 

29.4 

7.00 

50.0 

7.093 

143.3 

0.20 

5.000 

4.85 

33.3 

4.9 

50.0 

4.95 

100.0 

0.26 

3.846 

3.6 

15.63 

3.  7 

26.32 

3.  796 

76.69 

0.32 

3. 125 

3.05 

41.  7 

0.075 

62.5 

3.  1 

125.0 

0.38 

2.632 

2.5 

20.0 

2.582 

53.08 

2.6 

83.33 
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It  can  be  seen  in  Table  3 (or  from  X (t))  that  near  impact 
there  is  a large  increase  in  the  image  change  factor.  The  question 
to  be  answered  is,  what  value  of  X is  a limiting  one?  This 
illustrates  an  example  in  which  there  is  no  one  value  that  can  be 
specified  without  a full  prescription  of  the  system.  Rather  we 
would  postulate  several  conditions  and  then  enable  some  bound  to 
be  set. 

First,  the  initial  image  size  at  acquisition  is  determined 
by  the  range  and  optical  input  lens  focal  length.  Then,  the 
allowance  of  the  growth  of  the  image  size  expressed  through  the 
image  change  factor,  will  be  limited  by  the  transducer  format 
employed. 

Consider  the  focal  length.  The  maximum  image  change  factor 
will  occur,  for  a fixed  time  interval,  when  the  initial  image 
is  a minimum.  This,  in  turn,  will  be  established  by  the  resolu- 
tion limit  of  the  transducer.  A reasonable  value  for  the  MTF 
of  the  transducer  is  given  in  Ref.  4 as  20  cycles/mm  for  the  usually 
accepted  MTF  = 0.25  point.  Thus,  the  resolution  element  is 
0.05  mm  per  cycle.  If  we  also  accept  the  requirement  that  5 
cycles  are  required  across  the  target  for  detection,  the  minimum 
target  width  becomes  0.25  mm.  Using  Equation  (1),  the  minimum  ' 

focal  length  allowable  becomes  204  mm  (using  r"*"  because  it  leads  ’ 

to  the  smallest  image  at  the  centerline  acquisition  range  R^  = i 

1.525  km).  With  the  resolution-determined  minimum  image  size,  ! 

t 

the  transducer  will  be  filled  when  the  image  change  factor  j 

I 

becomes  144  and  228,  respectively,  for  36  x 24  mm  and  57  x 57 
mm  formats.  These  two  sizes  were  chosen  because  they  represent 
the  usable  portions  of  the  35  and  70  mm  film  formats. 
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Discussions  with  one  investigator  in  the  transducer  research 
area  indicate  that  these  two  sizes  are  technologically  acceptable 
and  that  sizes  larger  than  70  nun  may  lead  to  electric  field 
distribution  problems  and  thus,  require  further  investigation. 

We  could  establish  the  limit  on  the  image  change  factor  in 
another  way  by  prescribing  the  focal  length  first.  A typical 
value  would  be  305  mm,  a common  aerial  reconnaissance  camera, 
lens  focal  length.  With  this  lens,  the  limiting  image  change 
factors  are  26.08  and  42.47,  respectively,  for  the  35  and  70  mm 
formats.  In  these  cases,  the  centerline  slant  range  was  used. 
Several  computations  for  the  three  acquisition  ranges  are  given 
in  Table  4. 

TABLE  4 IMAGE  CHANGE  FACTORS  REQUIRED  TO  FILL 
I TRANSDUCER  FORMATS  (FL  = 305  mm) 


Initial 

Image  Change  Factor 

Acquisition 

to  Fill  Format 

Range  (km) 

35  mm 

70  mm 

= 2.985 

51.02 

83.09 

R = 1.525 
o 

26.08 

42.47 

r“  = 1.043 

17.84 

29.05 

Consider  now  the  two  lower  velocities,  214  and  305  m/sec 
for  the  same  acquisition  range,  the  minimum  range,  or  near  field,  ^ 

when  R^  = 1.525  km.  The  results  of  this  calculation  are  shown  || 

in  Figures  6 and  7 for  the  two  velocities,  respectively. 

Here  we  see  that  the  near  field  scene  image  change  is  a 

i significant  one  while  the  opposite  is  true  for  the  far  scene. 

I The  importance  of  this  difference  is  in  deriving  the  appropriate 

1 

1 i 
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matched  filters  for  the  scale  interval  to  be  covered  since 
the  initial  small  changes  for  the  far  scene  are  not  identical 
in  an  absolute  sense  to  those  for  the  near  scene. 

Although  the  upper  bound  on  the  image  change  factor  is  set 
by  the  maximum  image  that  can  be  placed  upon  the  transducer  with 
the  prescribed  optics,  the  difference  between  those  indicated 
in  Table  3 and  those  given  above  for  the  transducer  formats  does 
in  some  ways  establish  some  specification  for  the  zoom  lens 
requirements.  A chosen  format  established  the  maximum  allowed 
on  the  transducer;  the  desired  "last  look"  before  impact  establishes 
the  maximum  image  change  factor.  The  two  are  reconciled  by  the 
demagnification  of  the  zoom  lens  system.  Since  the  current  state 
of  the  art  for  these  lenses  shows  that  a lOX  magnification  can 
be  achieved,  we  can  expect  that  lOX  the  values  established  by 
transducer  limits  are  the  upper  limits  on  the  image  change  factor. 
Thus,  according  to  whichever  value  is  picked  from  the  above  dis- 
cussion, we  could  have  image  change  factors  of  261,  425,  1440, 
or  2280  depending  upon  the  format  and  focal  length  of  the  zoom 
lens.  We  see  from  Table  3 that  the  range  from  impact  potentially 
will  be  quite  small.  However,  before  the  benefits  of  this  optical 
advantage  can  be  taken,  one  must  also  reconcile  such  important 
influences  as  the  zoom  lens  speed  of  response  and  the  transducer 
update  capability.  Only  the  latter  is  considered  in  this  report. 

An  interesting  comparison  can  be  obtained  from  Eq.  (3). 

We  can  determine  the  variation  in  image  size  for  the  "near  field" 
scene  (i.e.,  for  R ) and  observe  the  simultaneous  change  in  the 
"far  scene"  (i.e.,  for  R^) . This  is  equivalent  to  the  projectile 
heading  toward  a close  target  and  yet,  having  another  target 
in  the  field  of  view  also  changing.  Of  course,  the  latter  target 
would  be  quickly  lost  from  the  scene  if  the  projectile  were  guided 
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to  the  closest  target, or  it  would  be  lost  from  the  field  of  view 
anyway  at  a later  time  as  the  areal  ground  coverage  reduces  as 
the  projectile  approaches  impact.  In  either  case,  we  can  observe 
how  the  relative  image  change  factors  for  both  scenes  are  altered, 
remembering  that  one  scene  is  lost  as  time  in  flight  progresses. 
These  comparisons  are  shown  in  Figures  6 and  7 for  the  velocities 
214  and  305  m/sec,  respectively. 

The  correlator  operates  by  storing  the  tactical  scene 
^ containing  the  target  and  then  playing  the  scene  through  the 

memory  bank,  one  element  of  which  contains  the  target  to  the  same 
' scale.  Each  time  the  transducer  recording  is  made,  it  is  antici- 

pated that  a matched  filter  has  been  recorded  to  the  same  scale 
and  can  be  comparable  to  the  input  image.  Considereing  then  the 
t foregoing  results,  we  can  see  that  there  are  significant  size 

changes  which  must  be  accounted  for.  Since  each  comparison  is 
done  in  a discrete  sequence,  we  must  determine  how  often  an 
individual  scene  is  recorded,  and  what  change  in  image  size  is 
allowable  before  the  individual  scene  is  recorded.  First,  we 
determine  how  many  images  are  required  in  the  memory. 


1 2 3 4 5 6 7 8 

TIME  (SEC) 


Fig.  6 Image  Changes  as  a Function  of  Time  for 
Target  Position  in  Scene  (V  ■ 0.21) 


5,  OPTICAL  MEMORY  BANK  REQUIREMENTS 


The  memory  density  requirements  depend  upon  the  image  size 
change  factor,  the  rate  at  which  update  can  be  achieved,  and  the 
parametric  sensitivity  that  is  employed  for  the  target  image. 

At  least  one  experience  (Ref.  6)  has  shown  that  the  (-3dB)  scale 
sensitivity  variation  is  ± 20  percent  (area)  for  the  M-60  in  an 
overhead  view.  The  linear  scale  variation  for  the  100  percent 
minimum  image  then  becomes  1.0  ± 0.1.  We  can  determine  the 
linear  (scale)  size  factor,  F,  for  any  variation,  that  prudence 
or  experimental  results  dictate,  from  the  fact  that  the  upper 
limit  for  one  matched  filter  is  the  lower  limit  for  the  adjacent 
one.  Thus, 


F 


1 -I-  g 
1 - a 


N 


(4) 


where  N « 0,  1,  2,  ...  is  the  number  of  filters  required  in 
order  to  achieve  a linear  size  factor  of  F when  a (symmetrical), 
-3dB  parametric  sensitivity  is  ± a.  The  coefficient  N must,  of 
course,  be  an  integer,  and  it  represents  the  whole  number  of 
matched  filters  needed  to  achieve  a linear  size  factor  of  F. 
Allowing  only  then  for  a scale  change  in  the  target  size,  we  can 
compute  the  number  of  filters  required  for  any  scale  (expressed 
relative  to  theinitial  size).  This  is  summarized  in  Table  5 for 
a ± 10  percent  scale  change,  -3dB  sensitivity.  We  can  carry  this 
analysis  further  and  consider  the  requirements  that  a scale 
change,  with  full  azimuthal  coverage  for  each  scale,  would  demand. 

The  criteria  for  the  allowance  for  orientation  variations 
are  also  derived  from  Ref.  6.  Although  the  data  are  based  upon  a 
strategic,  or  overhead  view,  and  the  present  application  requires 
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an  oblique  one,  the  data  to  be  used  are  believed  to  be  valid. 

This  is  based  upon  an  unpublished  investigation  which  showed 
that  the  average  sensitivities  of  a matched  filter  to  variations 
about  an  oblique  view  of  the  (tank)  target  are  comparable.  This 
means  that  head-on  views  were  more  sensitive  but  side  views  less 
sensitive.  Thus,  the  assumption  of  ± 8°  seems  to  be  a reasonable 
one. 

TABLE  5 NUMBER  OF  FILTERS  REQUIRED  FOR  + 10  PERCENT 
SCALE  CHANGE  SENSITIVITY 


If  we  now  assume  that  (1)  there  is  no  interaction  between  the 
parameters  scale  and  orientation,  and  (2)  the  orientation 
sensitivity  is  a constant  regardless  of  the  actual  scale,  then 
the  number  of  required  filters  can  be  obtained  from  Eq.  (4)  and 
a consideration  for  projectile  velocity.  Thus 


No.  of  Filters 
Required 

where  all  the  quantities  are  as  previously  given  and  W - 360/ 
(orientation  sensitivity  fullwidth)  taken  to  the  next  highest 
whole  number.  For  ± 8°,  W = 23. 

The  number  of  filters  required  for  a = ±0.1,  e . ± 8°  has 
been  computed,  and  is,  of  course,  dependent  only  upon  the 
parametric  sensitivities  and  the  largest  value  of  F allowed.  The 
results  for  the  case  cited  are  given  in  Table  6.  These  results 
are  a more  realistic  estimate  since  target  placement  and  size 
have  been  accounted  for. 

It  can  be  seen  that  many  positions  of  optical  memory  are 

required  for  representative  tolerances  in  the  parameters.  However, 

it  has  been  previously  shown  (Ref.  6)  that  recording  128  filter 

positions  on  an  area  of  40  x 40  mm  is  not  beyond  the  state  of  the 

art.  Through  such  techniques  as  overlay  and  interlace  this 

number  can  be  increased  to  500.  Also,  the  filter  capability 

cited  is  based  upon  a given  transform  lens  focal  length  (360  mm). 

Decreasing  this  would  increase  the  storage  capability  by  a factor 
2 

(360/new  FL)  for  the  same  area.  For  example,  we  could  obtain 
an  increase  in  storage  capacity  by  about  a factor  of  4 and  still 
be  above  the  storage  medium's  limiting  resolution.  Thus,  the 
seemingly  large  value  of  397  is  well  within  the  bounds  of  current 


log  F- 


+ 1 


(5) 
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TABLE  6 NUMBER  OF  FILTERS  REQUIRED  FOR  SCALE  CHANGE  OF 
+ 10  PERCENT  AND  ORIENTATION  SENSITIVITY  OF  ± 8 


technology.  It  should  be  pointed  out,  however,  that  the  task  of 
preparing  that  number  of  memory  positions  is  time  consuming  and 
tedious.  The  ameliorating  factor  is  that  usually  only  one  filter 
is  required  and  that  photographic  replication  (for  many  projectiles) 
might  be  possible. 

It  has  also  been  shown  that  the  degree  of  discrimination 
exhibited  by  a matched  filter  depends  upon  the  order  of  the  filter 
and  that  parametric  sensitivities  generally  associate  with  filter 
order.  Thus,  some  degree  of  parametric  latitude  should  be 
accounted  for.  This  has  been  done  by  using  Eq.  (4)  and  allowing 
a and  W to  assune  the  values  0.05,  0.10,  0.15,  0.20,  and  ± 8°, 

± 10°,  ± 12°,  respectively.  A few  selected  results  are  shown  in 
Table  7 for  the  size  factors  F = 26.1,  a value  obtained  from  the 
tran  ,ducer  -orii-aL;  filling  consideration  given  in  Table  4.  Similar 
results  could  be  computed  for  other  filling  factors. 

TABLE  7 NUMBER  OF  MATCHED  FILTER  MEMORY  ELEMENTS  REQUIRED 
FOR  VARIOUS  SCALE  AND  ORIENTATION  SENSITIVITIES 
(MAXIMUM  IMAGE  TO  FILL  TRANSDUCER) 


Scale 

Number  of 

Filters 

Required 

Parametric 

for  a - 

•3  db  Width  of 

Sensitivity 

+ 8° 

± 10° 

± 12° 

0.05 

772 

604 

504 

0. 10 

397 

311 

259 

0. 15 

271 

213 

177 

0.20 

208 

163 

136 

The  results  of  Table  7 clearly  indicate  that  OMFIC  memory 
requirements  are  a manageable  factor  in  the  system  concept.  This 
is  particularly  so  when  the  tactical  requirement  calls  for  zeroing 
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in  on  a class  of  vehicles  (e.g.,  tanks)  rather  than  a specific  one 
(e.g. , M-60  versus  M-48).  The  next  task  is  to  determine  how  often, 
and  how  fast,  the  transducer  must  be  recycled  in  order  to  meet 
the  requirements  dictated  by  the  foregoing  results. 
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6.  DETECTOR  RESPONSE- TIME  CONSTRAINTS 


It  is  apparent  now  that  one  of  the  unique  problems  associated 
with  the  guidance  of  projectiles  in  a terminal  homing  mission  is 
the  change  of  target  size  in  a fixed  focus  imaging  sensor.  If  we 
assume  that  the  final  detector  is  located  in  the  correlation  plane, 
i.e.,  it  detects  the  peak  correlation  between  the  target  being 
viewed  and  the  target  image  in  memory,  then  we  must  contend  with 
the  fact  that  the  peak  correlation  signal  will  vary  with  the  size 
or  scale  of  the  incoming  image.  One  way  to  avoid  the  impact  of 
the  problem  is  to  arrange  perfect  registration  of  the  correlation 
signal  from  all  the  target  sizes  from  the  different  positions  in 
the  matched  filter  memory  onto  a common  correlation  plane.  However, 
if  we  lost  registration  or  if  we  elect  for  signal- to- clutter 
considerations  to  have  separate  correlation  plane  detection  for 
each  memory  location  then  we  must  consider  the  time  variation  of 
an  individual  matched  filter  which  was  made  for  a fixed  range 
interval.  Since  we  must  sense  the  peak  signal,  this  would  impose 
a constraint  on  the  response  time  of  the  detector.  We  assume  in 
this  section  that  the  image  presented  to  t±ie  matched  filter  or 
correlator  has  not  been  smeared  by  any  other  components  such  as 
the  electro-optic  transducer  that  converts  the  input  incoherent 
image  into  a coherent  image.  This  is  the  case  when  the  response 
time  of  the  transducer  is  fast  enough;  the  burden  then  falls  on  the 
detection  plane.  A corresponding  analysis  of  the  effects  of  image 
scale  size  charge  is  given  later  in  this  report  to  address  the 
other  alternative,  i.e.,  where  the  transducer  becomes  the 
limiting  component. 

If  we  assume  that  the  target  is  imaged  by  a fixed  aperture 
lens,  and  the  projectile  is  moving  with  velocity  v then  the  linear 
image  size  can  be  expressed  as  (see  Eq.  (3)): 
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where  is  the  image  size  where  the  range  to  the  target  is  R 
and  the  elapsed  time  from  acquisition  t equals  zero.  We  note 
that  if  the  image  formed  by  the  sensor  at  time  t «*  0 is 
expressed  as  f(x,y),  where  x and  y are  coordinates  in  the  image 

2 

plane  and  f is  the  optical  amplitude  (the  intensity  in  watts/cm 
being  given  by  the  square  of  f ) , then  as  the  image  changes  size 
its  optical  amplitude  will  change.  If  we  assume  that  the  sensor 
has  a fixed  aperture  and  that  the  light  propagation  can  be 
computed  on  a spherical  wave  basis,  then  the  intensity  at  a 
given  point  in  the  image  plane  is  given  by 


(f/A)‘ 


(7) 


where  I is  the  intensity  of  the  equivalent  point  source,  F is 
s 

the  focal  length  of  the  system,  A is  the  linear  dimension  of  the 

aperture,  and  (f/#)  is  the  so  called  f-number  of  the  lens.  For 

a fixed  lens  we  see  that  the  point- by- point  image  intensity  is 

not  a function  of  range  to  the  target;  we  will  refer  to  this  as 

the  fixed  intensity  model.  Then  as  the  image  size  changes  the 

image  simply  expands  or  contracts  so  that  the  optical  amplitude 

of  the  image  is  given  by  f(ax,ay)  with  a = 1-vt/R.  If,  on  the 

other  hand,  we  assumed  that  the  light  propagation  from  the  target 

2 

is  not  spherical  with  a (1/R  ) intensity  fall  off  then  we  must 
consider  a variable  intensity  of  the  image.  To  bracket  our 
results  we  will  consider  the  case  where  the  total  light  flux  in 
the  image  remains  constant  as  the  image  changes  size.  This  then 
requires  that  the  image  intensity  decrease  as  the  image  grows 
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larger  to  keep  the  total  light  flux  constant.  This  case  yields 
our  optical  amplitude  given  by  a f(x,y)  so  that  the  integral 
of  the  optical  intensity 


2.2 
a f 


(ax,  ay) 


dxdy 


(r,s)  drds 


constant  (8) 


remains  the  same  as  the  image  changes  size  (or  a changes). 

The  matched  filtering  or  correlation  operations  on  the 
sensed  image  will  be  written  as 


C(^.Tl) 


// 


g(x,y)  h(x+a,  y+ii)  dx,dy 


(9) 


where  g is  the  sensed  image  (either  f(ax,ay)  or  a f(ax,ay) 
depending  upon  the  nature  of  the  imaging  system)  and  h (x,y) 
is  the  matched  filter.  We  will  consider  high  pass  matched  filters 
in  the  present  discussion  so  that  h(x,y)  is  a high  pass  replica 
of  the  given  image  f(x,y)  that  is  selected  at  some  nominal  scale, 
say  a » 1.  In  order  to  obtain  specific  results  we  will  take  a 
simple  rectangular  figure  as  the  prototype  f(x,y),  i.e. , 


f(x,y)  - ( 


l-i 


ly| 


L 

2 


(10) 


0 


elsewhere 


In  this  case  we  may  express  the  peak  value  of  the  correlation 
signal  C(^,ri)  as  (for  example  in  the  fixed  intensity  model) 


C(0,0) 


g(x,y)  h(x,y)  dx,  dy 
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■//' 


Cp  “ C(0,0)  I f(ax,ay)  f^p(x,y)  dx  dy 


f(ax,ay)  f(x,y)  dx  dy  - 


f(ax,ay)  f^p(x,y)  dx,  dy  (11) 


where  and  fj^p  refer  to  the  responses  of  the  high  pass  and  low 
pass  matched  filters,  respectively.  In  particular,  we  find  that 
that  if  a < 1 then  for  the  simple  rectangular  figure 


If 


il‘ 


f(ax,ay)  f(x,y)  dx  dy  =//  f (x,y)  dx  dy  - L L (12) 

X y 


and 


If 


# 


f(az,ay)  fj^p (x,y)dxdy  = ^^^2  I 


y [fj^p(x,y)]  dx  dy 


(13) 


*1 


' I 


we  then  have 


JJ^ E(ax, 


ay) 


(27r)- 


o)  L 

•2  SinJ^ 

CO 

X 

a 


CO  L 

2Sin 


CO 


CO  L , 

2 Sin 

2a 

CO 

a 


(15) 


d CO  d CO 

X y 


where  we  have  used  the  Fourier  transform  of  the  rectangle  and 
defined  B ,B  to  be  the  frequency  of  the  filter.  If  we  then  form 

y ^ 

a "figure  of  merit"  function  P given  by  the  ratio  of  the  square 
of  the  peak  correlation  function  at  various  scales  or  a' s to 
the  square  of  the  peak  correlation  function  at  a = 1 we  have, 
after  simplification 


in  the  fixed  image  intensity  model,  where 


B L 

X X 


Sin  — Sin  r 

a 


dr, 


1 

2 


B L 
V V 


Sin  ^ Sin  r 


dr 


B L 

X X 


o-  2 

Sin  r 


dr 


A 18) 


dr 


We  have  computed  a series  of  curves  to  indicate  the  expected 
variation  of  the  performance  function  P with  various  image  scales 


in  the  case  where  L 


2L  and  B 

X X 


B = B. 

y 


Figure  8 shows  the 


results  for  the  fixed  intensity  case  P^;  note  the  very  rapid 

B L BL 

changes  as  the  normalized  cutoff  frequency  x x x varies 

2*2 

from  t/4  to  ir  in  Figures  8(a)  and  8(b),  and  then  the  saturation 
effect  as  we  go  to  values  as  large  a 5v  in  Figure  8(c).  In  order 
to  illustrate  the  nature  of  the  effect  we  also  computed  the  value 
of  the  image  scale  change  that  would  produce  a performance 
function  value  of  1/2,  at  different  values  of  cutoff  frequency. 

In  order  to  ease  this  calculation  we  assumed  that  this  would  occur 
at  values  of  a near  1 so  that  we  could  write,  for  example,  a » 1 - 
and 


Sin  “ Sin  I + -y  Cos  I 

In  this  case  we  have,  continuing  the  first  order  analysis 

2 


(19) 


Pj-d 


key 


(20) 
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CONSTANT  INTENSITY 


where 


K 


dr  + A 


TT^ 


(21) 


1 - 


16 

2 


TT 


and  we  have  P^.  = 1/2  where 


vt 

R 


^3^4 


1- 


r2 


(22) 


These  results,  for  j^|<  0.1  are  plotted  in  Figure  9 together 
with  the  corresponding  results  where  we  consider  the  fixed  flux 
model  for  the  sensor.  We  see  the  saturation  effect  as  the 
normalized  cutoff  frequency  goes  beyond  2ir. 


As  an  example  of  the  impact  of  image  scale  change  on  the 

system  --  consider  the  case  where  V * 300  meters/sec  (Mach  1) 

and  R = 10  meters,  i.e.,  the  sensor  is  operating  at  a range  of 

10  meters  to  the  target;  then  if  we  utilized  a high  pass  matched 

BL 

filter  with  bandwidth  * t,  we  would  have 


f - 

for  a fixed  intensity  sensor  or 

t * 1. 7 msec 

However,  the  situation  is  highly  variable;  for  example,  if  the 

range  is  50  meters  and  the  velocity  is  150  meters/sec  (Mach  .5) 

we  have  T - 17  msec  to  read  the  correlation  signal  before  it 

falls  to  1/2.  That  is,  when  the  projectile  is  50  meters  from  the 

target,  it  takes  17  msec  for  the  output  correlation  signal  power 

to  drop  by  a factor  of  1/2  and  this  would  imply  the  need  to  have 

a detector  response  faster  than  17  msec. 
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7.  TRANSDUCER  UPDATING  REQUIREMENTS 


The  transducer  is  an  image  storing  device  consisting  of  an 
anisotropic  dielectric,  liquid  crystal  contained  between  layers 

■k 

of  light  sensitive  materials  and  transparent  electrodes.  When 
an  image  is  focused  upon  the  crystal  and  voltages  are  applied, 
the  image  is  recorded  in  the  alterations  of  the  crystal.  The 
application  of  other  appropriate  voltages  causes  erasure  and 
restoration  of  the  transducer.  The  elapsed  time  between  com- 
parable points  in  the  exposures  is  the  cycle  time,  and  this  time 
will  be  dictated  by  the  rate  at  which  one  wishes  to  record  and 
process  an  image. 

In  the  present  analysis  no  account  is  made  for  the  mechanics 
of  storage  or  erase,  or  other  aspects  of  the  transducer.  The 
concern  is  entirely  focused  upon  the  projectile  movement  and  the 
target's  image  upon  the  transducer.  Eventually  one  must  establish 
the  update  requirement  for  the  guided  projectile  system  on  the 
basis  of  the  guidance  system  and  the  specific  aerodynamics  of  the 
projectile.  Thus  the  number  of  guidance  updates  or  the  minimum 
distance  between  them  in  the  terminal  phase,  produces  the  require- 
ment as  the  transducer  update  or  cycle  time. 

The  present  discussion  is  not  specific  to  any  one  projectile 
and  thus,  we  cannot  make  a definitive  statement  as  to  the 
minimum  transducer  cycle  time  required.  However,  we  may 
arbitrarily  set  this  time  by  assuming  that  the  guidance  system 
will  update  each  time  the  optical  memory  address  is  changed  by 
an  amount  given  by  its  range  or  scale  sensitivity,  i.e.,  the  a' s 
introduced  previously.  Then  we  may  use  Eqs.  (3)  and  (4)  to 

The  discussion  applies  to  the  transducer  described  in  Refs.  4 
and  5;  a similar  transducer  is  described  in  Ref.  10  and  the 
reference  cited  therein.  Other  transducers  are  described  in 
Ref.  11. 
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establish  the  minimum  time  corresponding  to  the  last  update  at 
a fixed  scale  factor  change.  Equating  the  two  equations  leads  to 


where  as  before  V « v/R. 

Equation  (23)  has  been  computed  for  the  sequence  of  V-values 
for  the  scale  sensitivity  increments  a « 0,05,  0.10,  0.15  and 
0.20  for  two  cases.  The  first  is  for  times  quite  close  to 
impact  and  the  second,  for  the  case  where  the  transducer  is 
updating  at  close  to  television  frame  rates. 

The  results  for  the  first  case  are  summarized  in  Table  8, 
as  we  see  there  is  a requirement  for  update  rates  that  can  be  an 
order  of  magnitude  faster  than  TV  frame  rates.  The  columns  3-5 
i and  6-8  show,  respectively,  the  elapsed  time  and  the  duration  of 

I the  update  interval  for  the  case  where  N = 1,  and  the  case  where 

( 

! the  elapsed  time  is  close  to  the  impact  time  shown  in  column  9. 

Several  conclusions  emerge  from  Table  8.  Not  surprising, 
the  higher  V is,  the  shorter  the  update  interval  must  be  although 
it  can  also  be  seen  that  the  scale  parametric  sensitivity  can  be 
traded  off  for  velocity  ratio.  In  fact,  the  least  sensitive  filter 
operated  V - 0.38  requires  the  same  time  to  cycle  as  the  most 
' sensitive  filter  at  low  velocity.  This  comparison  is  significant 

j . ' when  we  consider  that  for  the  centerline  acquisition  range  « 

1.525,  V » 0.38  implies  a projectile  flying  at  Mach  1.7. 

For  the  more  realistic  ratios  V ■ 0.14,  0.20  the  transducer 
update  interval  average  is  about  7.5  msec.  This  is  about  an 


order  of  magnitude  shorter  than  current  technology  allows  for 
low  light  illumination  of  the  transducer.  For  bright  sunlight, 
recording  times  of  6 msec  have  been  demonstrated. 

The  results  of  the  second  case  -*  which  might  be  called  the 
current  technology  case  — are  summarized  in  Table  9.  The  column 
headings  are  similar  to  those  of  Table  8 except  that  an  additional 
column  has  been  added  to  show  the  range  from  the  target.  These 
data  are  significant  because  they  show  that  except  for  one  case 
(V  - 0.38,  a = 0.05)  the  transducer  performance  is  sufficiently 
satisfactory  to  update  to  a point  quite  close  to  the  target. 

Except  for  the  two  extreme  speeds  (V  « 0.32,  0.38)  this  distance 
is  d < 100  meters  and  therefore,  the  question  as  to  where  the  last 
guidance  update  is  given  becomes  pertinent  and  must  be  examined. 

Shown  in  column  7 is  the  range  "bin"  interval  which  is 
covered  by  the  memory  position  (-3dB  coverage)  corresponding  to 
the  parametric  sensitivity  increment,  a.  The  intervals,  of 
course,  vary  in  proportion  to  the  absolute  distance  from  impact. 
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TABLE  9 TRANSDUCER  UPDATE  REQUIREMENTS  - CURRENT  TECHNOLOGY 
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8.  IMAGE  ENLARGEMENT  SMEAR  - CONSTRAINTS 
IMPOSED  ON  TRANSDUCER  RESPONSE  TIME 


During  projectile  flight  the  duration  of  exposure  for 
recording  the  image  on  the  transducer  is  comparable  to  that 
for  photographic  exposures  in  aerial  reconnaissance  cameras.  It 
is  anticipated  however,  that  the  guided  projectile  may  not  have 
the  benefit  of  image  motion  compensation  (IMC)  now  a standard 
approach  in  cameras.  Also,  one  cannot  rely  upon  post  facto 
removal  of  motion  degradations  because  of  the  real  time  nature 
of  the  application.  Thus,  some  measure  of  image  motion  degrada- 
tion might  be  anticipated. 

In  aerial  cameras,  the  degradation  usually  takes  the  form  of 
a linear  image  motion,  defocus  or  vibration,  and  occasionally 
enlargement  smear.  In  the  OMFIC-GP  application  it  is  anticipated 
the  principal  cause  of  smearing  will  be  the  latter.  While  the 
presence  of  other  degradations  might  be  explored  for  an  implemented 
OMFIC-GP,  enlargement  smear  is  expected  to  be  present  and  the 
only  one  considered  here. 

If  a zoom  lens  system  is  employed  as  indicated  in  the 
schematic  of  Figure  3,  it  is  interesting  to  contemplate  the 
possibility  of  achieving  compensation  through  incremental  zoom- 
lens  control.  However,  this  implementation  problem  is  not 
exeimined  here  although  the  results  of  this  analysis  might  be 
useful  in  the  analysis  for  image  correction  through  a zoom  lens. 

We  can  use  Eq.  (2)  in  a modified  form.  If  we  consider  the 
instantaneous  image  rather  than  the  acquisition  image  size  on 
the  transducer,  we  can  then  interpret  the  image  as  the  "acquisition 
image"  when  the  shutter  is  open.  The  left  hand  side  of  Eq.  (2), 
designated  st,  represents  the  image  change  while  the  shutter  is 
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open.  Taking  differentials  of  £ and  t,  dividing  by  the  image 

size  when  the  shutter  is  open,  I , we  obtain 

s 


s 


100  = V 


(1  - Vt)‘ 


At  X 100 


(24) 


= 7o  Enlargement  Smear 


Equation  (24)  has  been  computed  and  the  results  are  shown  in 
Figure  10  for  two  velocity  ratios  (0.14  and  0.26)  and  for  exposure 
times  of  10  and  50  msec.  The  last  two  values  were  chosen  because 
they  bracket  anticipated  and  current  transducer  technology. 
Actually,  the  former  value  can  also  be  achieved  although  the 
illumination  requirements  call  for  bright  sunlight. 

It  is  not  surprising  that  the  shorter  exposure  results  in 
less  smear,  nor  that  near  impact  the  smear  is  greatest.  The 
latter  follows  from  the  results  described  in  Section  4 which 
showed  the  rapid  growth  in  image  size  near  impact. 

If  IMG  is  to  be  achieved  through  the  zoom  optical  system,  it 
is  clear  that  it  must  be  coordinated  well  with  the  system,  record/ 
read  cycle,  and  have  millisecond  reaction  times. 

In  the  remainder  of  this  section  we  will  evaluate  the  effect 
of  uncompensated  scale  change  on  the  correlation  signal  obtained 
with  a fixed  target  image  when  an  incoherent -to- coherent  optical 
transducer  is  present.  The  basic  effect  of  the  transducer  will 
be  to  smear  the  input  or  sensor  image  since  the  time  interval 
necessary  for  the  transducer  to  produce  sufficiently  intense 
images  might  be  so  long  that  we  would  have  to  consider  the  change 
of  image  size  during  the  exposure  time.  In  this  sense,  the 
transducer  is  acting  like  film  that  is  recording  an  image  whose 
size  changes  during  exposure  resulting  in  a blurred  or  smeared 
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RUNNING  TIME,  SECONDS 


Fig.  10  Enlargement  Smear  as  a Function  of 
Flight  Time  for  Two  Velocity  Ratios 


Image.  In  order  to  assess  the  effect  of  this  on  the  system  we 
will  compui.  the  effect  on  the  resultant  correlation  signal.  In 
Section  5 we  considered  the  effect  of  changing  image  size  on  the 
detection  of  the  correlation  signal  by  some  detector  with  limited 
response  time.  Depending  upon  the  type  of  components  to  be 
utilized  in  a particular  system  either  the  transducer  or  the 
detector  (or  possibly  both  if  their  response  times  are  equivalent) 
could  set  the  actual  system  limitation.  For  the  present  we  will 
however  treat  them  separately. 

If  we  assume  that  the  target  is  imaged  by  a ^ixed  aperture 
lens,  and  the  projectile  is  moving  with  velocity  v then  the  linear 
image  size  is  expressed  by  Eq.  (3).  If  the  image  is  formed  by 
a fixed  aperture  sensor  then  the  image  amplitude  can  be  written 
as 


f(ax,ay)  (25) 

2 2 
where  f is  the  image  intensity  in  watts /cm  ; the  model  assumed 

here  for  the  optical  system  is  that  of  spherical  wave  propagation 

with  1/R  intensity  falloff  and  the  image  intensity  I^  is 

determined  only  by  the  f/#  of  the  lens  (Eq.  (7)).  If,  on  the 

other  hand,  the  model  taken  for  the  optical  system  is  one  in  which 

the  image  has  a constant  flux,  i.e.. 


« constant 


then  the  optical  amplitude  is  given  by 


(26) 


af(ax,ay)  (27) 

and  the  integral  of  the  optical  density  given  by  Eq.  (8). 
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Now  In  order  to  model  the  transducer  we  will  assume  that  the 
resultant  transmission  of  the  transducer  after  exposure  is  given 
as  some  function  of  the  exposure,  i.e. , the  transmission  t is 
given  as  some  function  G(*)  of  the  total  optical  energy  density 
incident  on  it 


I 


g(x,y)  dt 


(28) 


where  T is  the  exposure  time  and  g(x,y)  is  the  optical  intensity 
falling  on  the  transducer  which  is  given  in  the  constant  flux 
sensor  model  by 


2 2 

g(x,y)  -of  (ox, ay) 
of  in  the  constant  intensity  model  by 


g(x,y)  - f (ax, ay) 


(29a) 


(29b) 


In  general,  we  could  have  an  arbitrary  functional  dependence  G(*) 
but  for  our  purposes  we  will  assume  that  the  image  intensity 
variations  are  small  and  we  adopt  a linearized  model  obtained  by 
expanding  the  expression  in  a power  series  in  y ; 


/’ 


thus  we  will  use 


t-wK  + K, 
o 1 


.y)  dt 


(30a) 


+ K,  / a 

o 1 J 


2 2 

f (ax, ay)  dt 


(fixed  flux)  (30b) 
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+ K,  / 
o 1 J 

o 


f^(ax,ay)  dt  (fixed  Intensity)  (30c)  A 


The  matched  filtering  or  correlation  operation  on  the  sensed 
image  after  it  is  recorded  on  the  transducer  will  be  given  as 


c(l,n) 


g(x,y)dt  h (x+0,(y+n)  dxdy  (31) 


where  g is  given  previously  and  h(x,y)  is  the  matched  filter's 
response.  We  will  consider  high  pass  matched  filters  in  the 
present  discussion  so  that  h(x,y)  is  a high  pass  replica  of  the 
intensity  of  the  given  image  f(x,y)  that  is  selected  at  the 
initial  range  R.  In  order  to  obtain  specific  results  we  will 
take  a simple  rectangular  figure  described  in  Eq.  (8).  We  also 
use  the  peak  correlation  signal  as  our  detection  signal;  for 
example,  in  the  fixed  intensity  case  we  will  use 


C = c(0,0) 


'//[/ 


y(x,y)  dt  fjjp(x,y)  dx  dy 


f^(ax,ay)  f^p(x,y)  dxdy  > dt 


' r f jf^(ax,ay)  |^f''(x,y)-f^p(x,y)J  dxdy  dt 
oL-'  (32) 
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In  particular.  If  we  take  the  situation  a < 1 we  have  for  the 
simple  rectangular  shape 


and 


// 


f^(ax,ay)  f^(x,y)  dxdy  - 


(33) 
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II 


f^(ax,ay)  f^  (x.y)  dxdy  - — ^ 
^ (27r)- 


O) 


o)  L\  / OD  L’ . 
2Sin  -hA  hsiXi  -ff' 


O) 


-B  "B  > 0)  L V / O)  L'  , 

>'  ’/2Sln-p\/2Sln-f^ 


U) 


Oi 


doa  do) 

X y 


(34) 


where  we  have  utilized  the  Fourier  transform  of  the  pulse  as 


and  defined  B and  B to  be  the  cutoff  frequency  of  the  filters. 

X y 

If  we  form  a "f igure-of-merit"  function  P given  as  the  ratio  of 
the  square  of  the  peak  correlation  signal  of  the  smeared  image 
at  various  exposure  times  T to  the  square  of  the  peak  correlation 
signal  of  the  constant  scale  (unsmeared)  image  at  the  same  exposure 
time  T then  we  have,  after  simplification 


(36) 
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in  the  fixed  image  intensity  model,  where  « “ and  the 
areas  previously  given  in  (Eq.  (16)). 

We  have  computed  a series  of  curves  to  indicate  the  expected 

variation  of  the  performance  function  P with  various  exposure 

times  in  the  case  where  L « 2 L and  B ■ B ■ B.  Figure  11 

y X X y 

shows  the  results  for  the  fixed  intensity  case  P ; note  the  very 

^ B L BL 

rapid  changes  as  the  normalized  cutoff  frequency  x x » x 

2 2 

varies  from  7r/4  to  tt  (Figs.  11a  and  11b)  and  the  saturation  as 

the  cutoff  frequency  extends  out  to  values  of  Sir  (Fig.  11c).  In 

order  to  illustrate  the  nature  of  the  effect  we  also  computed 

the  values  of  normalized  exposure  time  ^ that  would  produoe  a 

K 

performance  function  value  of  1/2  at  different  values  of  norn -lized 

BL 

cutoff  frequency  x . To  facilitate  this  calculation  we  assumed 

2 

that  this  would  occur  at  values  of  a near  1 so  that  we  could 

« 

write,  for  example,  a *=  1 - e and 


Sin 


2a 


Sin 


„ r 
2 + — Cos  2 


(37) 


In  this  case  we  have  (continuing  the  first  order  analysis) 


1 - . i I f dt 


1 _ k VI 
^ 2 R 


(38) 
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CONSTANT  INTENSITY 


Fig.  11b  Variation  of  Correlation  Signal  Due  to  Exposure  Time  (VT/R) 
Smearing  for  Constant  Image  Intensity  and  Various  Cutoff 
Frequencies  (BL/2).  2.5-^  in  t/4  steps 

(expanded  scale)) 


where 


k 

TT 


r sin  r , Tsin  r 

3/  ^ -T-^ 

Jo  •'O 


1 --^AA 
2 3 4 

TT 


In  particular,  Pj  “ where 


These  results,  forj'~j<  0.1  are  plotted  in  Figure  12  together 
with  the  corresponding  results  where  we  consider  the  fixed  flux 
model  for  the  input  image  to  the  transducer.  Again  we  see  the 
saturation  effect  as  the  normalized  cutoff  frequency  goes  beyond 
2t. 

As  an  example  of  the  impact  of  image  smearing  as  the  system 
we  consider  the  case  where  v « 300  meter /sec  (Mach  1),  R » 10 
meter,  and  we  utilize  a high  pass  filter  with  normalized  cutoff 
frequency  of  - tt.  From  Figure  12  we  have 


f-0.1 


so  that 


0.1  R 

t ■ “ *3.3  msec. 

V 

On  the  other  hand,  if  we  go  to  a range  of  R ■ 50  meters  and  a 
velocity  v = 150  m/sec  (Mach  0.5)  we  have 

t - 32  msec. 

That  is,  if  the  projectile  is  50  meters  from  the  target  and 
traveling  at  Mach  0.5,  then  the  image  smearing  on  the  transducer 
will  reduce  the  correlation  plane  signal  power  by  a factor  of  1/2 
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compared  to  the  signal  if  the  image  had  not  been  smeared;  this 
implies  the  need  to  have  a transducer  response  time  less  than 
32  msec,  to  avoid  this  penalty. 


I 


9.  OMFIC  IMPACT  UPON  A REPRESENTATIVE  SYSTEM  DESIGN 

Based  upon  the  preceding  analyses  one  might  consider  the 
"system"  aspects  of  the  OMFIC  configuration.  Specifically,  what 
performance  could  we  expect  in  a guided  projectile  (GP)  system? 

The  OMFIC  "guidance  system"  will  provide  a specific  target 
identification  (i.e.,  differentiation  of  an  M-60  from  a T-62) 
and  require  a relatively  modest  memory  bank  in  order  to  provide 
guidance  information  to  the  GP.  In  particular,  for  subsonic 
velocities  the  memory  bank  positional  requirements  will  be 
about  360. 

With  this  memory  bank  and  the  current  technology  transducer, 
the  GP  has  the  capability  of  providing  continued  guidance  for  the 
GP  from  an  initial  acquisition  range  of  1.5  km  to  within  approx- 
imately 50  meters  from  the  intended  target,  and  provide  a maximum 
range  uncertainty  of  ±5  meters  to  the  aerodynamic  guidance  system 
for  the  last  update. 

In  order  to  maintain  a correlation  signal  equal  or  greater 
than  -3dB  fromthe  peak  auto-correlation  value,  the  transducer 
must  "record"  an  Image  in  about  30  msec  or  less  (at  M^  ~ 0.5). 
This  is  somewhat  equivalent  to  having  the  transducer  operate  at  a 
standard  television  frame  rate,  a not  impossible  task  based  upon 
current  technology.  This  result  is  based  upon  an  uncompensated 
scale  change  effect  upon  the  correlation  signal  and  agrees  with 
the  situation  where  update  (transducer)  is  dictated  by  having 
a sequence  of  continuous  range  "bins."  However,  if  the  detector 
response  time  dictates  guidance  rates,  it  will  be  about  one-third 
faster  and  thus,  press  current  transducer  technology. 


Thus,  within  the  foregoing  limitations  the  OMFIC  system  is 
capable  of  providing  a GP  with  target  identification  and  projectile 
guidance.  In  some  projectile  (or  missile)  cases  it  is  understood 
that  the  aerodynamic  response  time  is  longer  so  that  the  OMFIC 
approach  becomes  more  viable. 
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